Magnetic excitations in multiferroic LuMnOs studied by inelastic neutron scattering 

H. J. Lewtas, A. T. BoothroydO M. Rotter, and D. Prabhakaran 
Oxford University Department of Physics, Clarendon Laboratory, Parks Road, Oxford, 0X1 3PU, UK 

H. MiiUer 
Institute for Solid State Physics, University of Technology Vienna, Wiedner Hauptstr 8-10, A-IO4O Wien, Austria 



o 

(N 

< 






M. D. Le 

Helmholtz-Zentrum Berlin fiir MateriaUen und Energie, 
Hahn-Meitner-Platz 1, D- 14109, Berlin, Germany 

B. Roessli and J. Gavilano 
Laboratory for Neutron Scattering, Paul Scherrer Insitute, CH-5232 Villigen PSI, Switzerland 

P. Bourges. 
Laboratoire Leon Brillouin, CEA-CNRS, CE Saclay, 91191 Gif sur Yvette, France 

(Dated: August 4, 2010) 

We present data on the magnetic and magneto-elastic coupling in the hexagonal multiferroic 
manganite LuMnOa from inelastic neutron scattering, magnetization and thermal expansion mea- 
surements. We measured the magnon dispersion along the main symmetry directions and used this 
data to determine the principal exchange parameters from a spin-wave model. An analysis of the 
magnetic anisotropy in terms of the crystal field acting on the Mn is presented. We compare the 
results for LuMnOa with data on other hexagonal RMnOa compounds. 

PACS numbers: 75.30.Ds, 75.25.Dk, 75.50.Ee, 75.85.+t 
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I. INTRODUCTION 

Multiferroic materials have been intensively studied in 
recent years following the discovery of compounds that 
display giant cross-coupling effects between magnetic and 
ferroelectric order parametersi"—. Particular interest has 
been aroused by the possibility of new magnetoelectric 
coupling mechanisms^^— and the potential for exploita- 
tion in technological applications^"—. One of the most in- 
vestigated families of multiferroics is the hexagonal man- 
ganites i?Mn03, which form with R = Sc, Y, Ho, Er, 
Tm, Yb and Lu. The magnetoelectric behavior found 
in this family is associated with frustrated antiferromag- 
netic interactions of Mn spins on a triangular lattice. The 
compounds formed with R = Sc, Y and Lu are attrac- 
tive for fundamental studies because they are not com- 
plicated by additional magnetic contributions from the 
R ions and because they form a family in which changes 
in structure and magnetoelectric behavior can be corre- 
lated in a systematic way. Here we focus on LuMnOs, 
and present neutron scattering measurements of the co- 
operative magnetic dynamics and measurements of the 
thermal expansion by dilatometry. The data provide 
quantitative information on the exchange interactions, 
magnetic anisotropy, and magnetostriction, all of which 
play a part in the magnetoelectric coupling. The results 
are compared with similar measurements on YMnOs. 

The crystal structure of the hexagonal manganites, 
which is described by the space group P63CTO, is built 
from corner-sharing MnOs bipyramids which form layers 
parallel to the ab plane separated by rare earth ions, as 



shown in Fig. [T] The Mn ions form a near-ideal triangu- 
lar lattice. The ferroelectric distortion, which occurs at a 
high temperature (Tc > 1000 K for LuMnOg, Ref.[l3), is 
caused by a tilting of the MnOs bipyramids and a buck- 
ling of the R plane, which together create a \/3 x \/3 su- 
perlattice distortion (trimerization) of the Mn ions and a 
ferroelectric moment along the c axisiSiii. The distortion 
shifts the Mn ions along the a axis away from the ideal 
X — 1/3 positioni^ — see Fig. [1] 

The magnetic properties of LuMnOa arise from the (al- 
most) triangular layers of Mn'^+ (3rf^) ions with S = 2. 
Neighboring spins are coupled by antiferromangetic ex- 
change interactions which are frustrated by the triangu- 
lar geometry, as evidenced by the large ratio of the Weiss 
to Neel temperatures |8/Tn| ^ 10 (Ref Il3l) . the anoma- 
lous magnetic entropy below T^^, and the reduction 
in the value of the ordered magnetic moment to about 
75% of the full spin-only moment for 5* = 2 (Ref. M)- 
The Mn spins in all the hexagonal i^MnOs compounds 
form a classical 120° structure within the triangular lay- 
ers (Fig. [T|>iSiiirJ^. The spins are confined by anisotropy 
to lie in the ab plane, and the large inter-layer separation 
decouples the layers electronically and makes the mag- 
netism quasi-two-dimensional. In the case of LuMnOs, 
magnetic ordering occurs below Tn = 88 K. 

Evidence for magnetoelectric coupling in i?Mn03 com- 
pounds is provided by anomalies at Tn in the di- 
electric constan t ^•^'^^ , lattice dynamics^Sr— , thermal 
conductivity^! and structural parameters'^"—. There are 
also interesting observations by optical second-harmonic 
generation which show a cross-correlation between fer- 
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FIG. 1: (Color online) (a) The hexagonal layered structure 
of LuMnOs. (b) Projection of the structure down the c axis 
showing the Mn sites. The large (red) diamond is the projec- 
tion of chemical/magnetic unit cell, and the filled and empty 
circles represent Mn at fractional heights z = and z — 1/2, 
respectively. The Mn trimerization has been exaggerated for 
emphasis. The in-plane (Ji, J2) and out-of-plane {J[, J'2) 
Mn-Mn exchange couplings used in the spin-wave model are 
shown. Arrows on the Mn atoms show the probable magnetic 
structure of LuMnOs based on neutron diffractionlfi and op- 
tical second harmonic generationi^. 



roelectric and magnetic domains due to the formation 
of magnetic domain walls below Tn which coincide with 
ferroelectric domain walls^""— . 

The precise microscopic mechanism of the magneto- 
electric coupling in i^MnOa has not been described yet. 
Careful structural measurements on Yi_a;Lua;Mn03 have 
shown that an isostructural transition takes place at Tn 
(Ref. 26) , which causes further displacements of the ions 
resulting in a small increase in the ferroelectric polariza- 
tion. It was therefore proposed that the magnetoelectric 
coupling is driven by a primary magnetoelastic coupling. 
The origin of the magnetoelastic coupling, however, re- 
mains unclear. One possibility is that the isostructural 
distortion may occur in order to relieve some energy as- 
sociated with magnetic frustration^. Another proposal 



is that the system might benefit energetically from the 
Dzyaloshinskii-Moryia interaction below Tn via a small 
c-axis displacement in the oxygen atoms that bond ad- 
jacent Mn atoms^-. This displacement would produce a 
small additional electric polarization along the c axis. 



II. EXPERIMENTAL 

Single crystals of LuMnOa were prepared by the op- 
tical floating-zone technique as follows. Polycrystalline 
LuMnOa was prepared by standard solid-state reaction 
from high purity (>99.999%) LU2O3 and Mn02. The 
polycrystalline powder was pressed into rods of diame- 
ter 8 mm and length 80 mm, and sintered at 1300°C for 
24 hours. Single crystals were grown in a four-mirror 
optical floating-zone furnace (Crystal Systems Inc.) at a 
scanning rate of 3-4 mm hr~^ with the feed and seed rods 
counter-rotating at SOrpm. The growth was performed 
in a flowing atmosphere of argon and oxygen in the ratio 
12:1. At each stage in the preparation the phase purity 
of the product was checked by powder X-ray diffraction. 

Unpolarized neutron scattering measurements were 
performed on a crystal of mass 1.9 g on the cold- neutron 
triple-axis spectrometer TASP at the SINQ facility (PSI, 
Switzerland) and on the thermal triple-axis 2T1 at LLB- 
Orphee (Saclay, France). At TASP, the crystal was 
mounted in an 'orange' helium cryostat and neutron spec- 
tra were recorded with a fixed final energy of 4.5 meV. 
The corresponding setup at 2T1 was with a closed-cycle 
refrigerator and a fixed final energy of 14.7 meV. Mea- 
surements were made with either a* and c* or a* and b* 
in the horizontal scattering plane, where a*, h* and c* are 
the axes of the hexagonal reciprocal lattice. On 2T1 some 
measurements were also made in the plane parallel to the 
a*h* plane but displaced by 0.5 reciprocal lattice units 
along the c* axis. The lattice parameters of LuMnOa 
referred to the space group PBacm are a = b — 6.05 A, 
c = 11. 4 A, and the inter-axis angles are a ^ (3 — 90°, 
7 = 120°. Hence, a* = b* = 47r/ (a^B), and c* = 27r/c. 

Magnetic measurements were performed with a super- 
conducting quantum interference device (SQUID) mag- 
netometer on a small piece of crystal cut from the same 
rod as the neutron crystal. The thermal expansion was 
measured on the same piece of crystal with a miniature 
capacitance dilatometer— 1^. 



III. RESULTS 

Figure [2] shows the magnetization of LuMnOa for a 
magnetic field of strength 1000 Oe applied parallel to 
the ab plane and along the c axis. Magnetic ordering 
is signalled by a sharp peak in the magnetization at 
Tn = 87.5 ± 0.5 K (Fig. [2] lower inset). This is confirmed 
by the appearance of magnetic Bragg peaks below Tn 
in neutron diffraction data (Fig. [5] upper inset). The 
magnetization exhibits a small anisotropy, being slightly 
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FIG. 2: (Color online) Zero-field-cooled magnetization mea- 
surements performed in a magnetic field of strength H = 
1000 Oe applied along the c-axis and a-axis. The lower in- 
set shows an expanded temperature range about the ordering 
temperature Tn = 87.5 ± 0.5 K. The upper inset displays the 
temperature dependence of the (101) magnetic Bragg peak 
amplitude. 



larger when the field is applied parallel to the ah plane 
{Xab) than along the c-axis. This easy-plane anisotropy 
is consistent with the observation that the moments lie 
in the plane in the ordered phase. The data follows a 
Curie- Weiss law at high temperatures (not shown) with 
a negative Weiss temperature, Q. From fits of 1/x vs T 
we obtain 8 = — 819±2K from Xab, and — 837±1K from 
Xc- These values are close to those reported previously 
from single crystals^^, but somewhat larger in magni- 
tude than obtained from powder samples'—. A clear cusp 
is seen at Tn in Xc whereas a broader peak is seen in 
Xab- It has been suggested that the cusp in Xc is caused 
by coupling between adjacent Mn layers, and the more 
rounded peak in Xab is due to frustration^'^. 

Examples of neutron scattering spectra from both in- 
struments are presented in Fig. |3l Figure [3l[a) shows 
energy scans recorded on TASP at the scattering vec- 
tors Q — (1,0,0) and (1,0,1.5), both of which con- 
tain two asymmetric peaks. Figure [3l^b) shows data at 
Q = (1.33,0.33,0.5) and (1.44,0.11,0.5) measured on 
2T1. Since the ordered moment on the Mn sites is rel- 
atively large (~ 3/ib) the scattering from magnons is 
expected to be much stronger than phonon scattering at 
these relative small scattering vectors. This, together 
with the resemblance of the spectra to previous measure- 
ments on YMnOa, gives us confidence that the main fea- 
tures in the spectra correspond to magnon excitations. 

To determine the magnon dispersion we fitted the 
peaks with Gaussian or Lorentzian functions (depending 
on the peak shape) on a linear background. The fitted 
peak positions have been collected together into a disper- 
sion curve which is plotted in Fig. 21 including data from 
both TASP and 2T1. Measurements on TASP extended 
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FIG. 3: (Color online) Neutron inelastic scattering from 
LuMnOs measured on (a) TASP at T = 5K and (b) 2T1 
at T = 13.5 K. The data are from constant-wavevector scans 
at the indicated positions in reciprocal space. The lines in (b) 
are fits to a lineshape comprising two gaussian functions. 



from below 1 meV up to about ISmeV, while measure- 
ments at 2T1 covered the range from about 5meV up 
to the energy of the highest modes. Measurements were 
performed in several different zones to find the maximum 
intensity and to check that the peaks were periodic in re- 
ciprocal space. 
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FIG. 4: (Color online) Spin- wave dispersion of LuMnOa. 
Solid circles are the peak centers obtained from fits to scans 
such as those shown in Fig. [S] Solid lines are calculated 
from the model described in the text with parameters Ji = 
-4.09(2) meV, J2 = -1.54(5) meV, J2 = +0.019(2) meV, 
,][ — 0, and D = — 0.48meV. The inset is a sketch of the 
{h, k, 0) plane in reciprocal space showing the path FMKF. 



Figure [S] shows diffuse scattering measurements in the 
vicinity of the point (1,0,0) in reciprocal space at tem- 
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FIG. 5: (Color online) Diffuse neutron scattering around 
the forbidden nuclear reflection (100). The inset shows the 
temperature variation in the diffuse scattering intensity at 
(1,0,0). 



peratures close to Tn- The (100) reflection has zero nu- 
clear structure factor, and we observed zero magnetic 
intensity here at low temperatures. The scans in Fig. [5] 
reveal strong difi^use scattering at temperatures close to 
Tn = 87.5 ± 0.5 with maximum intensity at Tn itself, 
as shown in the inset. The diffuse scattering is highly 
anisotropic, being very broad in the (0, 0, /) direction but 
much sharper in the {h, 0, 0) direction. 
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FIG. 6: (Color online) Thermal Expansion of LuMnOs mea- 
sured along the hexagonal a- and c directions. The inset shows 
the thermal expansion coefficient a as determined from the 
temperature derivative of the strain dL/L. 



Figure[6]shows measurements of the thermal expansion 
of LuMnOs parallel to the a and c axes. A magnetoe- 
lastic anomaly is clearly visible at the Neel temperature, 
both in the strain dL/L and in its temperature derivative 
a = L~^dL/dT. On cooling through Tn the magnetoe- 
lastic strain expands the c axis and shrinks the hexagonal 
plane. As a check, we also measured the thermal expan- 



sion in the hexagonal plane in the direction normal to a. 
The data resemble the behaviour of da/a to within the 
experimental error. Above Tn the thermal expansion is 
highly anisotropic. The c-axis strain is almost tempera- 
ture independent. 



IV. ANALYSIS AND DISCUSSION 

We first review the magnetic structure of the .RMnOa 
compoundsi^^— , with particular reference to LuMnOs. 
The six Mn sites in the unit cell (Fig. [Ij form two near- 
equilateral triangles, one in the z = layer and the other 
in the z — 1/2 layer. The spins on these triangles lie in 
the basal plane and order in 120° structures. Symmetry 
constrains the relation between the z = and z = 1/2 
layers to two possibilities, conventionally labeled a and 
/3. In the a structure the spin on the Mn at (x, 0, 0) {x w 
1/3) is parallel to that at (1 — x, 0, 1/2), whereas in the 
/3 structure the spins on these two sites point in opposite 
directions. Each spin makes an angle ip to the unit cell 
axis on which it lies, and the magnetic structure factors 
depend on ip and on the stacking relation (a or j3). For 
X — 1/3, the magnetic structures occur in homometric 
pairs, such that the magnetic diffraction intensities for 
the configuration {a, ip) are identical with those from (/3, 
V'±90°). 

In the case of LuMnOa, Katsufuji et al. concluded 
from neutron powder diffraction measurements^'' that 
the low temperature structure is one of two possibilities, 
either (a, ip = ±90°) or (/3, V == 0° or 180°). These 
two structures transform respectively like the r4 and 
r2 irreducible representations of the space group PG^cm 
(Ref. KLSr ) . Only this particular homometric pair have a 
completely absent (100) magnetic reflection. In our single 
crystal measurements we also found very little intensity 
at the (100) reflection (but relatively strong intensity for 
the (101) reflection — sec the inset in Fig. [2]) at the low- 
est temperature (T = 2K), in agreement with Katsufuji 
et al. The homometric pairs can in principle be distin- 
guished by optical second harmonic generation (SHG). 
Using this method, Fiebig et al. found that their sample 
of LuMnOa was a two-phase mixture of a structures with 
■0 = 0° and an unspecified other ip valuei^. Since our low- 
temperature data, as well as that of Ref. [l^ conclusively 
rule out any a structure which does not have tp ~ ±90° 
it is difficult to see how to reconcile the diffraction and 
SHG results. For the purpose of modeling the spin wave 
spectrum we will assume the (a, ip = ±90°) structure, as 
shown in Fig. [T] The magnetic spectrum of this and its 
homometric partner are not distinguishable at the level 
of precision of our data. 

We calculated the spin wave spectrum from the spin 
Hamiltonian 
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with two in-plane near-neighbor interactions ( Ji and J2) 



and two inter-plane interactions (J{ and J 2) defined as 
shown in Fig. [T] The first summation in ([T]) is over pairs 
of spins with each pair counted once so that the J con- 
stants are per spin pair. The second term models the 
out-of-plane anisotropy with a single-ion anisotropy pa- 
rameter D. We neglect the small in-plane anisotropy 
since the in-plane magnon gap was too small to measure 
in our experiment. 

Analytic expressions have been given previously for the 
spin- wave energies derived from spin Hamiltonians simi- 
lar to (H)), Refs. 36-39. These expressions have been ob- 
tained via the usual transformation of the Hamiltonian 
into a quadratic form of boson normal-mode operators 
in the linear approximation. Here we use an alterna- 
tive method based on dynamical matrix diagonalisation 
(DMD) as outlined in previous work^, which is imple- 
mented in the McPhase software package^. This for- 
mulation employs the random phase approximation to 
calculate the magnon cross sections in addition to the 
dispersion relations. 

There are six Mn spins per unit cell, which gives rise 
to a total of six spin- wave modes for each wavevector. As 
the interlayer coupling is small, the in-plane dispersion 
relations appear as three branches each containing two 
nearly-degenerate modes. The degeneracy of the lowest 
two modes is lifted close to the F point and along FZ 
by the effect of the J{ and J 2 interactions, while the 
upper four modes are almost degenerate along FZ. This 
degeneracy precludes the possibility to fit accurate values 
for J[ and J2 independently, and so we chose to fix J( — 
and to vary J2 under the constraint that Jj > to 
maintain the stability of the a structure. 

A least-squares fitting procedure returned the follow- 
ing values for the model parameters: Ji = —4.09(2) meV, 
J2 = -1.54(5) meV, J'2 = -f 0.019(2) meV {J[ = 0), and 
D — — 0.48meV. The calculated dispersion relations from 
the model with these parameters are shown in Fig.|4l The 
agreement is seen to be very good, and the parameters are 
well constrained by the data. For example, the 6.5 meV 
gap to the upper mode at F is sensitive to the single-ion 
anisotropy, and J2 controls the dispersion in the out-of- 
plane direction. The splitting of the magnon peaks in the 
vicinity of the K point, seen in Fig. [31 is sensitive to the 
difference between Ji and J2- Qualitative agreement be- 
tween the measured and calculated magnon cross sections 
gave us further support for the obtained parameters. 

As a check, we calculated the bulk magnetization us- 
ing the best-fit exchange and anisotropy parameters. The 
T/v predicted by the mean-field model is about 2.5 times 
larger than the observed Tjv, presumably as a conse- 
quence of frustration. The calculated susceptibility has a 
small easy-plane anisotropy consistent with the measured 
susceptibility. Fig. [5J The magnetization of the sample 
as a function of applied field (not shown) is linear and 
almost identical in the a and c directions. Increasing the 
single-ion anisotropy in the model creates a step in the 
a-axis magnetization which is not observed. This adds to 
the evidence that the single-ion anisotropy is very small 



compared to the exchange interactions. 

The exchange parameters obtained here show that the 
dominant magnetic interaction is the in-plane antiferro- 
magnetic superexchange via the ~120° Mn-O-Mn path. 
The inter-layer superexchange is two orders of magnitude 
weaker, confirming that the magnetism in LuMnOa is 
highly two-dimensional. It is interesting to compare the 
magnetic spectrum of LuMnOa investigated here with 
those obtained from similar measurements on YMnOa 
(Refs.milliiii) andHoMnOg (Ref.lssI). Qualitatively, 
the spectra of the three compounds look very similar, but 
the overall band width of the LuMnOa spectrum is about 
30% larger than that of YMnOa and HoMnOa (21meV 
compared with 16meV). Consistent with this, the fitted 
exchange parameters for LuMnOa are found to be sys- 
tematically larger than those of YMnOa and HoMnOa- 
This accounts for the difference in the antiferromagnetic 
ordering temperatures of these compounds: Tn ~ 88 K 
(LuMnOa) compared with Tn ~ 72 K (YMnOa) and 
Tn ~ 75 K (HoMnOa) and the larger Weiss tempera- 
ture of LuMnOa (6 « -850 K) compared with YMnOa 
(e«-700K)i3. 

The stronger magnetic interactions in LuMnOa fits 
with the systematic trend in the ionic radii and the lat- 
tice parameters^, i.e. the smaller the ionic radius the 
smaller the unit cell and the stronger the exchange inter- 
actions. The single-ion anisotropy parameter D is also 
found to be larger for LuMnOa {D = — 0.48meV) than 
for YMnOa {D ^ -0.28 to -0.33 meV) and HoMnOa 
{D — — 0.38meV). This could be another consequence of 
the small differences in the structural parameters of these 
compounds. 

The character of the diffuse scattering from LuMnOa 
close to Tn (Fig.[S]) strongly resembles that observed from 
YMnOa.— ^ The appearance of scattering which is very 
broad along c but relatively sharp in the plane indicates 
that the inter-layer correlations are very weak, consistent 
with the small J[ and Jj and two-dimensional nature of 
the magnetic system. The strong enhancement in the 
diffuse scattering intensity around (1,0,0) was also ob- 
served in powder diffaction measurements on LuMnOa J^ 
These showed that the diffuse peak persists up to at least 
~ 3Tn^ which was interpreted as evidence for strong 
geometric frustration. 

Careful powder diffraction measurements^ have 
shown that the magnetically-induced ferroelectricity in 
i?MnOa is associated with an isostructural transition 
involving an additional rotation of the MnOs bipyra- 
mids, and that the increase in ferroelectric polarization 
below Tn scales with the square of the ordered mo- 
ment. The thermal expansion of LuMnOa reported here 
(Fig. [S]) reveals a striking magnetoelastic anomaly at 
Tn, consistent with the diffraction data of Lee et al. 
(Ref. [23) who argued that magnetoelastic coupling (ex- 
change striction) is the primary source of the magneto- 
electric coupling^^. One might expect, therefore, that 
the magnetoelectric effect would scale with the strength 
of the exchange interactions and hence be greater in 



LuMnOa than in YMnOs. Support for this idea is 
provided by the magnetically-induced polarization cal- 
culated from the measured atomic displacements below 
Tn, which indeed appears to be systematically larger for 
LuMnOa than for YMnOa (Ref. 261. However, given the 
large experimental uncertainties in the values of the small 
atomic displacements this evidence should be considered 
tentative. 

Another notable feature of the thermal expansion is 
how small the c-axis thermal expansion is relative to the 
ab plane expansion above Tn (see Fig. ^. This effect, 
which is observed both for LuMnOs and YMnOs, does 
not correlate with the compressibility of these materi- 
als, which is similar along the c direction and in the ab 
plane^'^ . The relatively isotropic compressibility suggests 
that the anomalous c-axis thermal expansion is not due 
to a straightforward anharmonicity in the interatomic po- 
tentials along the c axis, and it would be interesting to 
find out what is responsible for it. 

Before concluding, we investigate the origin of the 
small magnetic anisotropy of LuMnOa, which is repre- 
sented in the Hamiltonian ^ by the phenomenological 
D{S^Y term. For reference, we performed a point-charge 
calculation of the the crystal field at the Mn sites assum- 
ing Mn to be in the Mn^"'" state with d^ configuration and 
using the structural parameters reported by Katsufuji et 
al. in Ref. [lOl. We included only the five nearest oxy- 
gen neighbors of Mn in the Mn05 bipyramid, as shown 
in Fig. [71 The ground-state S = 2 manifold is split by 
the crystal field via the spin-orbit interaction. The cal- 
culated splitting, due to the crystal field alone, is also 
shown in Fig. [T] This model of the crystal field predicts 
that the ordered magnetic moment points along the nor- 
mal to the local mirror plane, as shown in Fig. [I] and 
with the inclusion of the exchange field from the neigh- 
boring Mn ions predicts an anisotropy gap in the magnon 
spectrum of ^lOmeV. This size of gap is in clear disagree- 
ment with the observed spin-wave modes, which have a 
gap of less than 1 meV at the zone centre, Fig. |31 The 
anisotropy gap may be reduced in the model if the local 
symmetry is increased to Cat,, i.e. by lessening the degree 
of Mn trimerization and tilting of the MnOs bipyramid. 
We conclude, therefore, that the single-ion anisotropy is 
controlled by a tiny distortion of the ideal MnOs bipyra- 
mid and that the anisotropy is much smaller than that 
predicted by a simple point-charge model. 



V. CONCLUSIONS 

We have measured the magnon dispersion in LuMnOa 
and achieved a very good description of the spectrum 
using the Heisenberg Hamiltonian, Eq. ([TJ. We also ob- 
served a striking magnetoelastic coupling at Tn in the 
thermal expansion. The results are qualitatively simi- 



lar to those previously obtained on the sister compound 
YMnOa- The bandwidth of the one-magnon spectrum 
of LuMnOa is about 30% larger than that of YMnOa, 
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FIG. 7: Left: Calculation of the thermally-averaged charge 
density in the MnOs bipyramid of LuMnOa at T = 4K. The 
3d charge is represented (in blue) by a surface of constant 
charge density obtained from a calculation of the crystal field 
acting on Mn^+ assuming point charges of — 2|e| on each of 
the five nearest oxygen neighbors (shown in red). Right: Low- 
lying energy levels of Mn"^"*" split by the point-charge crystal 
field via spin-orbit coupling. The symmetry of each level 
is labeled according to the irreducible representation of the 
point group Ca, which describes the local symmetry around 
the Mn site. 



and the difference between the two nearest-neighbour 
in-plane exchange constants Ji and J2 is greater for 
LuMnOa than for YMnOa. As the magnetic interac- 
tions are stronger in LuMnOa than in YMnOa we ex- 
pect the magnetically-induced ferroelectric polarization 
to be greater in LuMnOa. The available diffraction data 
provides tentative support for this. 
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